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Abstract. We study the random acceleration model, which is perhaps one of 
the simplest, yet nontrivial, non-Markov stochastic processes, and is key to many 
applications. For this non-Markov process, we present exact analytical results for the 
probability density p{tm\T) of the time tm at which the process reaches its maximum, 
within a fixed time interval [0,T]. We study two different boundary conditions, which 
correspond to the process representing respectively (i) the integral of a Brownian bridge 
and (ii) the integral of a free Brownian motion. Our analytical results are also verified 
by numerical simulations. 
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1. Introduction 

Consider a general stochastic process x{t), starting from x(0) = 0, over a fixed time 
interval [0,T]. Let tm denote the time at which the process achieves its maximum value 
Xm during the interval [0, T] (see Fig. 1). Clearly, tm is a random variable that fluctuates 
from one realization of the process to another. Our main goal is to investigate the 
probability density function (pdf) p{tm\T) of the stochastic times tm, given the interval 
length T and the underlying stochastic process x(t). 

This question naturally rises in a variety of fields. For instance, in the context of 
the queueing theory, the stochastic process x{t) may represent the length of a queue 
at time t, and one would like to know tm, i.e., the time at which the queue length is 
maximum over a given time span. In the context of finance, x{t) may represent the 
price of a stock, and a trader is evidently interested in tm, i.e., the time at which the 
stock price is at its highest. 

Perhaps, one of the simplest examples is provided by the underlying stochastic 
process x{t) being an ordinary Brownian motion 

|^,W. (1) 

where 7]{t) is a Gaussian white noise with zero mean {rjlt)) = 0, and a delta correlator 
{rj{t)r]{t')) = 2D6{t — t'). In this case, the probabiliy density p{tm\T) is independent of 
D and is well known [1] 

p{tm\T) = ^ :; 0<tm<T. (2) 

ytm{T — tm) 

This curve is symmetric around the mid-point tm = T/2, has a '?7 shape' with minimum 
at tm = T/2, and diverges at the two end points tm = and tm = T. The cumulative 



distribution Prob(tm < x\T) = |arcsin ytm/T is known as one of Levy's celebrated 
'arcsine laws' pj. 

For a Brownian bridge, i.e., a Brownian motion constrained to return to at T, 
i.e., x{T) = 0, the corresponding pdf is known to be uniform: p{tm\T) = 1/T for 
0<tm<T ^. 
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Motivated principally by the applications in the queueing theory and finance, the 
pdi p{tm\T) has been computed explicitly for a variety of other 'constrained' Brownian 
motions by using suitably adapted path integral methods. These examples include 
Brownian excursion ^ , Brownian meander [3] , refiected Brownian bridge [3] , Brownian 
motion till its first-passage time and Brownian motion with a drift [5]. Some of these 
exact results were then reobtained via a functional renormalization group method [6]. 
In Ref. |6], the authors computed p{tm\T) also for Bessel processes and for continuous- 
time random walks. Curiously, the pdf p{tm\T) also appeared recently as an important 
input in the calculation of the area enclosed by the convex hull of a planar Brownian 
motion of duration T [7[ E], which displays an interesting application in ecology. In 
addition, we have recently shown that p{tm\T) also appears in the computation of the 
disorder-averaged equilibrium distribution of a particle moving in a random self-affine 
potential [9]. 

The results for p{tm\T) mentioned above have been obtained for Brownian motion 
and its variants, which are all Markov processes. The purpose of this paper is to go 
beyond Markov processes and present an exact result of p{tm\T) for a non-Markov 
process. The non-Markov process that we study here is the well-known random 
acceleration process, which evolves via 

x{t) = r/(t), (3) 

rjit) being a Gaussian white noise with {ri{t)) = 0, as before, and {'q{t)rj{t')) = 2D6(t—t'). 
For simplicity, we will choose D = 1 subsequently. The process starts at x(0) =0 with 
f (0) = x(t) = 0, and evolves until the time T. We denote by xj = x(T) and Vf = v(T) 
the final position and velocity of the process, respectively (see Fig. 1). The random 
acceleration process in the single variable x{t) is non-Markovian, whereas the position- 
velocity process {x, v} is Markovian and does not depend on the past history. 

The random acceleration process, being among the simplest non-Markovian 
processes, has been intensely studied both in Physics and Mathematics literature. 
In Physics, for instance, it appears in the continuum description of the equilibrium 
Boltzmann weight of a semifiexible polymer chain with non-zero bending energy |10j . 
It also describes the steady state profile of a (1 + l)-dimensional Gaussian interface [TT] 
with dynamical exponent z = 4, the continuum version of the Golubovic-Bruinsma- 
Das Sarma-Tamborenea model [T2]. This process also appears in the description of 
the statistical properties of the Burgers equation with Brownian initial velocity |13j . 
The first-passage and a variety of other related properties of the random acceleration 
model are highly nontrivial and have been studied extensively over the last few 
decades [lin5l[T6l[l0l[I7l[l9l[IIl[l8l[20l[2T]. 

Thus, in addition to being relevant in many applications, the random acceleration 
model represents a simple, yet nontrivial, non-Markov process where one can try to 
compute observables of physical interest. Recent studies have concerned the distribution 
of extreme observables associated with this process, notably the global maximum Xm 
itself over the interval [0, T] [TOl [22l |23]. In this paper, we focus instead on the time 
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Figure 2. A realization of a random acceleration process reaching its maximum 
at tm- In this case, the maximum occurs at the end of the time interval [0,r]. 



at which the global maximum Xm occurs in [0,T]. Actually, we show in this paper that 
p(tm\T) can be computed exphcitly. 

Let us first summarize our main results. Since the only time scale in the problem 
is T, it is evident that p(tm\T), normalized to unity over < tm < T, has the scaling 
form 

p{tm\T) = ^pC-^), (4) 

where the scaling function p{z), defined over < 2; < 1, satisfies the normalization 
condition: p{z)dz = 1. We will consider two different boundary conditions, detailed 
below, for which we are able to compute p{z) explicitly. 



1.1. Integral of a Brownian Bridge 

When the final velocity vanishes, Vf = 0, the process can be interpreted as the integral 
of a Brownian Bridge. In this case, we will show that 

^^"^ ^ Am[zii-z)f'' 

which is evidently symmetric around the mid-point z = 1/2 and diverges at the two end- 
points as 2;"^/^ and (1 — 2;)^^^^, respectively. It is also useful to consider the cumulative 
distribution P{z) = Jq p{z') dz', which reads 

^W-^B.fiiV (6) 



r2(l/4) n4'4. 

Here Bz{p,q) = Jq x^"^ (1 — xY^^ dx is the incomplete Beta function [21]. A plot of 



P{z) is shown in Fig. 4, where it is also compared to direct numerical solutions, to an 
excellent agreement. 



1.2. Integral of a free Brownian Motion 

When Vf is arbitrary, the process can be interpreted as the integral of a free Brownian 
motion. In this case, we obtain the following exact result for the normalized pdf 

p{z) = C6{z-l)+ ^^^z-'/' (1 - z)-'l\ (7) 
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where the constant C has the exact value 

C =1- ^ = 0.387628 ... (8) 

The density is thus asymmetric around the mid-point tm = T/2 (i.e., z = 1/2), and the 
maximum may either occur at some time strictly shorter than T, namely tm <T (i.e., 
z < 1), or with a finite nonvanishing probability C = 0.387628.. at the end point of the 
interval t^ = T (or equivalently z = 1). The representative trajectories for these two 
situations are shown respectively in Fig. 1 and Fig. 2. In other words, roughly 38.67% 
of all trajectories, starting at a;(0) = and f (0) = 0, achieve their maximum only at the 
end of the interval [0,T]. The corresponding cumulative distribution P{z) = Jq p{z') dz' 
is given by 

P„).C<,(.-l) + ii^ft(i5). (9) 

where 9{z — 1) is zero for z < 1 and is equal to 1 for 2; = 1, i.e, P{z) exhibits a 
discontinuous jump at z = 1 from 1 — C = y^3/8 = 0.612372 ... to 1. A plot of P{z) 
is provided in Fig. 5, where it is also compared to direct simulation results: we find an 
excellent agreement between analytical and numerical results. 

The paper is organized as follows. In Section 2, we outline the main ideas used 
for the exact derivation of p{tm\T) via path decomposition techniques. In Subsection 
2.3, we also compare our analytical predictions with the results obtained via Monte 
Carlo simulations. A concluding Section 3 presents a summary and raises some open 
questions. The details of the calculations are left to four Appendices, as they involve 
rather cumbersome multiple integrations of Airy functions. 

2. Calculations 

The basic ingredient in our computation is the propagator Z+(xi, f 1; Xq, fo, i), i.e., the 
probability that the process x(t), starting at a;o > with velocity fo, reaches the point xi 
with velocity vi at time t, without ever crossing the origin x = during this time. The 
Laplace transform of the propagator has been explicitly computed by Burkhardt and is 



recalled in Appendix A We show that the location of the maximum can be expressed 
in terms of the propagator Z+. Then, we make use of the known results about Z+ in 
Laplace space to derive a closed form expression for p{tm\T). 

The velocity v{t) of a random acceleration process is a Brownian motion, thus 
v{t) is continuous everywhere, which implies that x{t) is continuous and differentiable 
everywhere. An important consequence is that the global maximum Xm can lie either 
inside the interval (0,T), with a velocity f m = (Fig. 1), or at the boundary tm = T, 
with a velocity Vm = Vf > (Fig. 2). 

With reference to Fig. 1, the global maximum of the process starting at 
{xq = vq = 0} and ending at {xf,Vf} lies inside the interval (0,T). It is useful to 
decompose the total path in two portions: a first path from {x(0) = 0,f(0) = 0} to 
{xm,Vm = 0} in a time tm, and a second path from {xm,Vm = 0} to {xf,Vf} in a 
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Figure 3. a) Left. Decomposition of a realization of the random acceleration process 
into two paths, from to tm (I) and from tm to T (II). b) Right. Illustration of the 
change of variables x — Xm — x and t = t„i — t for (I) and x — Xm — x and t ^ t — 
for (II). 



time T — tm (see Fig. 3a). Remark that Xm represents the maximum of both paths. 
Denote by 'Plxm.Vm = 0;x(0) = v{Q) = 0,tm] the probabihty of the first path and by 
P[xf,Vf;Xm,Vm = 0, T — tm] the probabihty of the second path. With respect to the 
variables {x, v}, the process is Markovian, and the probabihty of the entire path is given 
by the product 

P[Xm,Vm = 0;x(0) = f(0) = 0,tm] X P[x f , V f] Xm, Vm = 0,T - tm]- (10) 

We introduce the change of variables x = Xm — x, t = tm — t (which implies v = v) 
for the first path (see Fig. 3b). Therefore, we can rewrite 

F[xm, Vm = 0; x(0) = v{0) = 0, tm] = Z+{xm, 0; 0, 0, tm). (11) 

As for the second path, the change of variables (see Fig. 

3b) gives 

F[Xf,Vf;Xm,Vm = 0, T - t J = Z+{Xf = Xm- Xf,-Vf;0,0,T -tm). (12) 

When tm = T, as illustrated in Fig. 2, an analogous argument, with a change 
of variables and t = T — t (which implies v = v), allows rewriting the 

probability of such a path as Z+{xm, 0; 0, vj, T). 

2. 1 . Integral of a Brownian Bridge 

In this case, Vf = 0, and the location of the maximum lies inside the interval, as in Fig. 
1. It follows that p(tm\T) can be obtained by integrating the two paths of Eq. ffTOj) over 
Xm and Xf, with f/ = 

dxm J dxfF[xm, 0; 0, 0, tm]F[xf, 0; 0, T - t J = 

roo roo 

/ dxmZ+{xm,0;0,0,tm) dx f Z +{x f , 0] 0, , T - tm) = 
Jo Jo 

l2{0,tm)l2iO,T-tm). (13) 
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We define 

l2{e,t) = / dxZ+{x,0;e,0,t), (14) 
Jo 

which denotes the probability that the process, starting at the initial position e > at 
t = 0, with initial velocity Vq = 0, remains positive up to time t, with a vanishing final 
velocity (the final position being arbitrary). 

The integral J2(e, t) can be computed exactly and it turns out that, in the e — ?■ 
limit, this integral vanishes, as apparent from Eq. flC.ll|) . This is not surprising: indeed, 
if the particle starts at the origin with vanishing velocity, it can not survive up to a 
finite time t, since it will cross the boundary almost immediately. This is due to the 
continuous nature of the Brownian velocity. Therefore, if one puts e = straightaway, 
the rhs of Eq. f|T3l) vanishes. The reason for this is clear: the Ihs of Eq. f|T3|) represents 
the probability that the maximum lies in [tm,tm + dtm], i.e., p{tm\t)dtm, and hence is 
proportional to the small time increment dtm- Setting e — ?• essentially implies dtm = 0, 
which therefore gives on the rhs. To extract the nonzero probability density p{tm\t), 
one therefore needs to keep a finite e (therefore a finite dtm) in /2(e, t) on the rhs of 
Eq. f|T3|) . Then, finally one uses the following limiting procedure that gives a finite 
answer 

We note that this type of regularization procedure has been adopted before in computing 
the distribution of the area under a Brownian excursion [25] and, more recently, in the 
context of computing p{tm\T) for several constrained Brownian motions, including the 
Brownian excursion [3j. 

2.2. Integral of a free Brownian Motion 

When the final velocity f/ of the particle is arbitrary, there is a finite probability that 
the global maximum occurs at the end of the observation time window [0,T]. This 
means that the probability density p(tm\T) includes a delta function C6(tm — T) (with 
a nonzero weight C ). This is in addition to a non-delta function part with a nonzero 
support (and a total weight (1 — C)) over the full interval tm G [0, T]. 

Let us first compute the delta-function component in p{tm\T) at tm = T. This 
contribution comes from paths that start at the origin with zero velocity and end up at 
Xm (with Xm being the maximum) in a small time window tm G [T — 5, T], where 5 — ?■ 
and Xm is free to take any positive value. Following the notations and the change of 
variables discussed earlier, one can write the net probability of such paths as 

lim/ p{tm\T)dtm= / dXmdvfZ+{xm,0;0,Vf,T) = Ii{T). (16) 
s^oJt-s Jo Jo 



We compute this integral h{T) explicitly in Appendix B and it turns out to be a 
constant indepedent of T, h{T) = C = 1 — ySjs. Thus, it follows that the delta- 
function component of the probability density is p{tm\T) = C6{tm — T). 
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Figure 4. Simulation results for the cumulative distribution P{z) = p{z') dz' 
(circles) as compared to the analytical formula in Eq. (|6]) (solid line), for the integral 
of a Brownian Bridge. 



We next focus on the non-delta function part, where the maximum occurs at some 
time tm well inside the interval [0,T]. The calculation is performed along the lines of 
the integral of a Brownian Bridge. The integral of the two paths of Eq. ( ITOl) now writes 

/ dXmZ+{xm,0;0,0,tm) / dv f / dxfZ+{xf,-Vf;0,0,T -tm) 

Jo J-oo Jo 

= 12(0,0/3(0, T-tJ, (17) 

where 

/+00 roo 
dv / dxZ+{x,v;e,0,t). (18) 
-00 Jo 

The integral /3(e,t) is computed in Eq. (ID. 71) and, for reasons explained already, it 
vanishes as e — )■ 0. Using the regularization mentioned before with xq = e and taking 
into account the full normalization, /^^ p{tm\T) dtm = 1, we can then write 

pitm\T) = C5{tm -T) + {l-C) hm -^2(e,t )/3(e,r-U ^^g^ 

^ ' ^ ^ ^ ^ '^^^jUtmh{e,tm)h{e,T-tm) ^ ' 

This thus gives the result of Eq. ([7]) . 
2.3. Numerical simulations 

To verify our main theoretical predictions in the two cases, namely, (i) the integral of a 
Brownian bridge in Eq. ^ and (ii) the integral of a free Brownian motion in Eq. dH]), 
we have also performed Monte Carlo simulations of the two processes. In both cases, 
we simulated 10"^ realizations, with T = 1 and an integration step 10~^. In Fig. 4 
and 5, for the two cases (i) and (ii) respectively, the theoretical and numerical results 
are compared. The circles represent the simulation results and the solid lines represent 
the analytical formulae, respectively in Eq. (jH]) (Fig. 4) and Eq. ([9]) (Fig. 5). The 
agreement between simulations and analytical curves is excellent. 
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Figure 5. Simulation results for the cumulative distribution P{z) (circles) as 
compared to the analytical formula in Eq. ([9]) (solid line), for the integral of a free 
Brownian motion. 

3. Summary and conclusions 

In this paper we liave studied two non-Markov stochastic processes: (i) the integral of 
a Brownian bridge up to time T and {ii) the integral of a Brownian motion up to T - 
i.e., the so called random acceleration processes. For both processes, we have computed 
exactly the probability density p{tm\T) of the time tm at which the process achieves its 
maximum when observed over a time window [0,T]. In the former case, the density 
p{tm\T) is symmetric around the midpoint tm = T/2, with power law divergences 
at the end points tm = and t^ = T. In contrast, in the latter case, p{tm\T) is 
asymmetric around the midpoint tm = T/2 and, in addition, has an unusual delta- 
function contribution at the end point tm = T. The exact results are given respectively 
in Eqs. ([5]) and ([7]). Our analytical findings are in excellent agreement with numerical 
Monte Carlo simulations (see Figs. 4 and 5). 

Note that, even though our final results look rather simple, their derivations, 
starting from the basic propagator, are nontrivial and involve rather delicate 
mathematical manipulations of multiple integrals, which we have tried to present 
systematically in a stepwise fashion in the Appendices. Hopefully, some of these details 
might be also useful in other problems involving the integrals of a Brownian motion. 

The exact result for p{tm\T) that we obtain in this paper is directly relevant 
for a particle diffusing in a one dimensional random potential V{x). In a finite box 
of size T, the particle finally thermalizes with an equilibrium Boltzmann density, 
Pcq{x) = e~^'^^^y Z , where (5 is the inverse temperature and Z = e~^^^^^ dx is the 
partition function. In the low temperature limit, by taking the average over disorder, 
one can show |9] that Peq{x) is precisely equal to the probability density p{tm = x\T) 
that the potential V{x), viewed as a stochastic process in space (where the space plays 
the role of 'time'), has its minimum (or equivalently its maximum [26]) at tm = x. For 
example, for the well known Sinai model, where V{x) itself is a free Brownian trajectory, 
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Peq{x) = l/liry x{T — x)], the same as the arcsine law. Our finding for p{tm\T) then 
generahzes the resuh for Pcqix) to the case where the potential V{x) represents the 
trajectory of a random acceleration process. 

Let us end with a curious observation. For the free Brownian motion, Levy's 
arcsine law appears in the distribution of two different observables: (a) in p{tm\T), 
i.e., the location of the maximum, and (b) in p{toccup\T), where toccup = lo i^i''')] 
represents the occupation time, i.e., the time the process spends on the positive half 
axis within the interval [0, T]. One may then ask if these two observables share the same 
distribution also for the random acceleration process. While we are able to compute 
p{tm\T) explicitly in this case, the computation oip{tocc\ip\T) seems harder in the random 
acceleration case. Our numerical evidence shows that, unlike in the Brownian case, the 
two distributions are different in the random acceleration model. Therefore, computing 
the occupation time distribution p(toccup|7') for the random acceleration model remains 
a challenging open problem. 



Appendix A. Propagator with an absorbing boundary 

In the absence of boundaries, the free propagator for the random acceleration process 
reads ITOl 



Zo{x,V]Xo,Vo,t) 



73 r 3 



In presence of the absorbing boundary at a; = 0, the propagator reads [TO] 

Z+{x,v]XQ,VQ,t) = Zo{x,v;xo,Vo,t) + Zi{x,v; Xo,Vo,t), 
where the Laplace transform of Zi{x,v;xo,Vo,t), i.e., 

Z,(x,.;.„,..,.) = /Ae-%(...;.,..„,i), 
has a rather complicated expression [TO] 

Zi{x, v; xo, vo, s) = -— / / dGdFAi{s/G^/^ - vG^^'^) 

2tt Jo Jo 

exp {- [Fxo + Gx+ |s3/2 + G^^)] } 



{x - Xo - Vt) (x - Xo - Vot) + — {v - Vq 

o 



(A.2) 
(A.3) 



Ai(s/F2/3 + t;oFV3)- 



{FGyi^F + G) ' 

Ai(x) being the Airy function [27]. In subsequent calculations, we will also use the 
following amazing identity [161 [10] 



27r 



°° dG 



F + G 



^-1/6 



exp 



2 /s^ ^ 
3 



Ai 



/ s 



VG2/3 



ir-i/6Ai 



JP2/3 



Appendix B. The first integral 

In this Appendix we compute the first integral 

Ji(T) = / / dvdx [Zo{x,0;0,v,T) + Zi{x,0;0,v,T)] 
Jo Jo 



(B.l) 
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where the propagator Zq is given in Eq. (lA.ip , and the Laplace transform of Zi is given 
in Eq. f lA.4|) . We show below that /i(T) is actually a constant, independent of T, and 
is given by an amazingly simple expression 



h{T) = C 



3 

1 - \ -. 



(B.2) 



Performing these integrals in closed form requires the use of a number of identities from 
Ref. [21]. However, to use the available identities, we need to first express the integral 
in a suitable form. To summarize, the derivation is not straightforward and requires 
quite a few mathematical manipulations and tricks. To make it easy, we lay out the 
main steps used in this derivation. 

Step 1: Upon substituting the exact Zq from Eq. (lA.ip . the first integral in Eq. (IB. ID 
can be performed in a straightforward manner and gives 

5 



poo roo ^ 

/ / dv dx Zq{x^O;0,v,T) = —. 
Jo Jo 12 



(B.3) 



Step 2: To perform the second integral, we first take its Laplace transform with respect 
to T, substitute the result from Eq. (]A.4ll . and then easily perform the integration over 
X. This yields 

poo poo poo 

J{s)= / / dvdxdTe-'^Zi{x,0]0,v,T) -- 
Jo Jo Jo 



27T Jo 



oo roo roo 



(FG) 



-1/6 



G{F + G) 



exp 



2 
3 



,3/2 



+ 



;3/2^ 



Ai 



Q2/3 



Ai 



/ s 



dv dG dF 
(B.4) 



Step 3: Next, we split the term 1/{G{F + G)} in the integrand on the rhs of Eq. (]B.4p 
into two parts, namely, 

(B.5, 



G{F + G) F{F + G) FG' 

and thus split the integral J(s) into two parts, i.e., J(s) = Ji{s) + J2(s). This allow 
dealing each contribution separately, so that a number of identities can be used, as 
explained below 

Step 4: Let us first consider the first term, 

1 



Ms) . 

{FGy^ 

F{F + G) 



oo roo POO 







exp 





2 /s 

3 [ 



dv dG dF 

3/2 3/2 \ 



Ai 



Q2/3 



Ai 



+ vF^/A. (B.6) 



F G 

Now, the integral over G can be performed explicitly using the identity (lA.Sp . Next, for 
the integral over v, we make a change of variables: w = F^^^v + sF~'^^^. Substituting 
back, we get 



Ms) 



dF 



Ai 



F2/3 



Ai(w) dw. 



(B.7) 
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Then, making a further change of variables, sF = z, yields a simpler expression 

3 roo 



^ POO POO 

Ms) = — dzAi{z) / Ai{w)dw. 
2s Jo Jz 

The rhs can be computed using integration by parts, so to give 



Ms) 



3 

4^ 



Ai{w)dw 



n2 



Using Ai{w)dw = 1/3 [27], finally gives the rather simple exact expression 
- ik- 



(B.8) 



(B.9) 



(B.IO) 



Step 5: We now turn to the second contribution J2(s), which is given by 

dv dG dF 



1 roo roo poo 

J2{S) = -- / / 

ZTT JO JO Jo 

exp 



GiF + G) 



2 /s3/2 s3/2^ 
3 



Ai 



Q2/3 



Ai 



/ s 



+ vF^/A . (B.ll) 



F G 

Now, the integral over G can be separated from the integrals over F and v, and we write 

1 



■Ms) 



2tt 



J2l{s) J22{s) 



where J2i{s) involves only the integral over G 

2 



J2^(s) 



POO 

/ dGG-^/^exp 
Jo 



G 



Ai 



Q2/3 



(B.12) 



(B.13) 



and J22{s) involves the integration over F and v 

2 



J22S) 



POO 

/ dF F-^l^ exp 
Jo 



/ Ai{F^/^v + sF-'^'^)dvm.U) 
Jo 



Step 6: We now perform the integral J2i{s) in Eq. (lB.13p . We first make a change of 
variables from G to y = 2s^/^/3G. Next, we use the identity 



Ai 



2 / 



TT 



1 /3\i/3 



V3 \2 



(B.15) 



where Kyiy) is the modified Bessel function of index u [21]. Substituting back in 
Eq. flRT3|) . we get 

1 1 /""^ 

The integral over y can be explicitly performed [21], and gives 7r-\/27r. Thus, we get a 
very simple expression 



J21 fsl 



^1/4' 



(B.17) 



Step 7: Now we turn to the computation of J22{s) in Eq. (IB. 141) . Making the change 
of variables from v to w = F^/^v + sF~'^^^, we get 

2 .00 
3\~F^j JsF-y^ 



POO 

Ms) = / exp 

JO 



Ai(w)(itf . 



(B.18) 
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Next, we make a change of variables from F to y = 2s^/^/3-F, so to rewrite the integral 

as 



1 



33/4 V 2 Jo 



dyy 



-1/2 ^-s/ 



(3j//2)2/3 



Ai{w)dw. 



Substituting the expressions for J22{s) and J2i(s) in Eq. flB.12p finally yields 



J2(S] 




An 



where the constant Aq is given by the integral 



An 



dyy 



^1/2 g-y 



(3y/2)V 



Ai{w)dw. 



(B.19) 



(B.20) 



(B.21) 



Step 8: We have now to evaluate the constant Aq in Eq. (]B.21|) . As a first step, let us 
again use the identity in Eq. (IB.lSp . so to express the inner integral in Eq. (lB.2ip as 



1 f°° 

Ai{w)dw = — = / ^1/3(^)^2;. 



, ..w^,.,^.. (B.22) 

This follows by making a change of variables w = {3z/2)'^^^ on the Ihs of Eq. flB.22p and 
then using the identity in Eq. fIB.lSp . Thus, we get 

1 



roo roo 

Ao = ^ dyy-'/'e-y Ky,{z)dz. 
7rv3 "'0 Jy 



(B.23) 



Next, we perform integration by parts and use the incomplete Gamma function defined 
by r[a,?/] = dx x""'^ e'^ . This gives 

1 



An 



vr / Ki/s{z)dz 



K,/,{y)T[l/2,y] dy 



Using the result [21] Kii^{z)dz = n/VS, we get 
3 ttVS Vo 



A. 



Ki/3(y)r[l/2,y] dy 



3 TTV^ 

Step 9: Concerning finally the integral Ai in Eq. (IB. 251) . namely, 

POO 

Ai= / Ki/3(y)r[i/2,y]f/y, 



we first write 



r[i/2,y] = -^e-^-ir[-i/2,y]. 



(B.24) 



(B.25) 



(B.26) 



(B.27) 



which can be easily proved by integration by parts. Substituting this result in Eq. (]B.26|) . 
the first term can be easily computed, and gives7rv^27r. We then get 

1 



Ai = 7rV27r 



2 Jo 



K,/,iy)T[-l/2,y] dy. 



To perform the second integral, we first use the integral representation 



T[-l/2,y] = —y-'/'e-y 



TT 



y + t 



dt. 



(B.28) 



(B.29) 
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Substituting this result in Eq. (IB.28P gives 

VTT JO Jo y + t 

Next, we use the identity [2 



y + t yt 



(B.30) 



(B.31) 



so to perform the inner integral over y in Eq. (lB.30p . The integral over t then becomes 
simple and using once again Ki/^{t)dt = vr/y^ yields an exact expression for Ai, 
namely, 



Substituting this result in Eq. (IB.25P gives the final expression for Aq 
Ao=\l 



TT. 



(B.32) 



(B.33) 



Step 10: Using the expression for Aq in Eq. (]B.20|) we compute J2{s), which, combined 
with the result for Ji{s) in Eq. (IB.lOp . allows evaluating J{s) = Ji{s) + J2{s) 



J{S) = J,{S) + J2{S) 



7 
12 



(B.34) 



This shows that the inverse Laplace transform of J(s) is just the constant prefactor of 
1/s in Eq. ( lB.34p . Combining this result with Eq. (1B.3P finally yields our main result 



h{T) = C=l 



(B.35) 



Appendix C. The second integral 

roo 

h{e,t)= / da;[Zo(a;,0;e,0,t) + Zi(a;,0;e,0,t)]. 

The first term can be easily integrated, and expanding in small e yields 

1 ^3 



dxZo{x, 0; e, 0, t) = —= 

Ay/TCt 2ttV 



e + ... 



(C.l) 



(C.2) 



We next take the Laplace transform of the second term with respect to t: 
Zi{x, 0; e, 0, s) = Zi{x, 0; e, 0, t) e~^^ dt and then integrate over x using the expression 
in Eq. ( lA.4p . To extract the leading e dependence, it turns out to be useful to split 
1 + 1. This gives 



dxZi{x, 0; e, 0, s) 



/ dF(e 
Jo ^ 



1 $(F,s),(C.3) 
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where 

1 r'^ dG 



2 fs 

exp 



Ai ( tIh) Ai f-l^ ) . (C.4) 



2 ^ 

The integration over F in the first term of Eq. (]C.3P can be carried out by resorting to 
the identity (1A.5I) . so to get 

m s)dF = - r T^Ai^ (-^) = --!=. (C.5) 



Jo 6-4/3 VC^/sy 4v/i' 

This contribution, after performing the inverse Laplace transform, cancels the one 
coming from Eq. (]C.2|) (at the leading order in the small e expansion). 

We address now the second term in Eq. (1C.3I) in the small e expansion. By means 
of the change of variables eF = z, the leading order of the second term reads 

dF (e-^^ - l) ^F, s) = ^ ^ ' B,B2is), (C.6) 

^ ^ ZTT 



where 



and 



Bi = I dz^-^ = -r(-l/6) (C.7) 
B.(s) = l ^exp --— A,(^). (C.8) 



Setting y = s'^^'^/G, we get 

1 dy _ ( 2 \ ^ 2/3\ 

/o ?/ 

Using Ai (0) = 1/ (3^/31 (2/3)), we finally get 

l-,F{e~---,)HF,s) = '^^g,. (C.IO) 

Then, by inverting the Laplace transform, we obtain the leading order contribution 

Appendix D. The third integral 

Remark that the integral 

/■ + 00 l-OO 

I'i{e,t) = I dv dx [Zo{x,v; e,0,t) + Zi{x,v; e,0,t)] (D.l) 







represents the survival probability at time t of the process started at the origin with 
vanishing velocity. The Laplace transform of this quantity, with respect to t, has been 
computed in [10], and reads 



isie, s) = -- I dF^{F, s) - I dF{e-'' - 1)^(F, s), (D.2) 
s Jo Jo 



oo pca 
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where 



Ai s/F 



72/3 



FV3 



1 2 ,s'^ '2 , 
2' 3 F ' 



40F 



(D.3) 



For the second term, we can set s/F"^/^ — y, so to make the dependence on s exphcit, 
and get after some algebra 



/•oo Q i-oo 

/ dF^{F,s)^- dyki{y) 
Jo Zs Jo 



1 + 



1 

s 



(D.4) 



This contribution cancels the ^ of the first term. The third term can be computed by 
resorting to the substitution z — eF, namely 



f dFie-^^-lMF, s) = eV3Ai(0) f (e^ - l) 
Then, by noting that, for small x, r{—l/2,x) we get 



-n/l 2 es^/^ ' 
1 _|_ ^ 2' 3 2 ^ 

40F 



(D.5) 



dF{e-'' -l)^{F,s) 




871^'^^Ky^Jo zV^ 



By finally performing the inverse Laplace transform, we obtain 

25/6r(_4/3) ei/6 



32/37rr(3/4) iV4 



+ 0(eV3). 



25/6r(-4/3) e'/"^ 
327% ^ 



(D.7) 



,(D.6) 



References 



[1 

[2: 

[3: 
[4; 
[5; 
[6; 
[7: 
[s: 

[9 

[lo; 
[11 
[12: 



Levy, P.: Sur Certains Processus Stochastiques Homogcnes. Comp. Math. 7, 283 (1939) 

Feller, W.: An Introduction to Probability Theory and its Applications. New York, Wiley (1968) 

Majumdar, S. N., Randon- Fur ling, J., Kearney, M. J., Yor, M.: On The Time To Reach Maximum 

For A Variety Of Constrained Brownian Motions. J. Phys. A: Math. Theor. 41, 365005 (2008) 
Randon-Furling, J., Majumdar. S. N.: Distribution of the time at which the deviation of a Brownian 

motion is maximum before its first-passage time. J. Stat. Mech. P10008 (2007) 
Majumdar, S. N., Bouchaud, J. -P.: Optimal Time to Sell a Stock in Black-Scholes Model. 

Quantitative Finance 8, 753 (2008) 
Schehr, G., Le Doussal, P.: Extreme value statistics from the Real Space Renormalization Group: 

Brownian Motion, Bessel Processes and Continuous Time Random Walks. arXiv:0910:4913 
Randon-Furling, J., Majumdar, S. N., Comtet, A.: Convex Hull of N planar Brownian Motions: 

Application to Ecology. Phys. Rev. Lett. 103, 140602 (2009) 
Majumdar, S. N., Comtet, A., Randon-Furling, J.: Random Convex Hulls and Extreme Value 

Statistics. J. Stat. Phys., in press, arXiv: 0912:0631 
Majumdar, S. N., Rosso, A., Zoia, A.: Hitting Probability for Anomalous Diffusion Processes. 

Phys. Rev. Lett., in press, arXiv:0911:3815 
Burkhardt, T. W.: Semiflexible Polymer in the Half Plane and Statistics of the Integral of a 

Brownian Curve. J. Phys. A: Math. Gen. 26, L1157 (1993) 
Majumdar, S. N., Bray, A. J.: Spatial Persistence of Fluctuating Interfaces. Phys. Rev. Lett. 86, 

3700 (2001) 

Golubovic, L., Bruinsma, R.: Surface diffusion and fluctuations of growing interfaces. Phys. Rev. 
Lett. 66, 321 (1991); Das Sarma, S., Tamborenea, P.: A new universality class for kinetic growth: 
One-dimensional molecular-beam epitaxy. Phys. Rev. Lett. 66 325 (1991) 



17 



[13] Valageas, P.: Statistical properties of the Burgers equation with Brownian initial velocity. J. Stat. 
Phys. 134, 589 (2009) 

[14] McKean, H. P.: A winding problem for a resonator driven by a white noise. J. Math. Kyoto Univ. 
2 227 (1963) 

[15] Goldman, M.: On the first passage of the integrated Wiener process. Ann. Math. Stat. 42 2150 
(1971) 

[16] Marshall, T. W., Watson, E. J.: A drop of ink falls from my pen... it comes to earth, I know not 

when. J. Phys. A: Math. Gen. 18, 3531 (1985) 
[17] Sinai, Ya. G.: Distribution of some functionals of the integral of a random walk. Theor. Math. 

Phys. 90, 219 (1992) 

[18] Lachal, A.: Last passage time for integrated Brownian motion. Stoch. Proc. Appl. 49, 57 (1994) 

[19] Majumdar, S. N., Bray, A. J.: Persistence With Partial Survival. Phys. Rev. Lett. 81, 1142 (1998) 
[20] De Smedt, G., Godreche, G., Luck, J. M.: Partial survival and inelastic collapse for a randomly 

accelerated particle. Europhys. Lett. 53, 438 (2001) 
[21] Burkhardt, T. W.: Dynamics of Absorption of a Randomly Accelerated Particle. J. Phys. A: Math. 

Gen. 33, L429 (2000) 

[22] Gyoergyi. G., Moloney, N. R., Ozogany, K., Racz, Z.: Maximal height statistics for 1/f" signals. 

Phys. Rev. E 75, 021123 (2007) 
[23] Burkhardt, T. W.: First-Passage and Extreme- Value Statistics of a Particle Subject to a Constant 

Force Plus a Random Force. J. Stat. Phys. 133, 217 (2008) 
[24] Gradshteyn, I. S., Ryzhik, I. M.: Tables of Integrals, Series, and Products. Academic, New York 

(1980) 

[25] Majumdar, S. N., Comtet, A.: Airy Distribution Function: From the Area under a Brownian 

Excursion to the Maximal Height of Fluctuating Interfaces. J. Stat. Phys. 119, 777 (2005) 
[26] Note that for a symmetric process, the distribution of the location of its minimum coincides with 

that of its maximum. 

[27] Abramowitz, A., Stegun, I. A.: Handbook of Mathematical Functions. Dover, New York (1965) 



